function fun_createFigure(F, A, B, freq, I_row, I_col, ...
        noise_intensity, noise_mean_num, ...
        phase_shift, phase_shift_tilts, filepath)
    
    % 生成均匀相移的矩阵
    Figure = cell(F,1);
    Delta = cell(F,1);
    for f = 1 : F
        Delta{f,1} = phase_shift(f)*ones(I_row,I_col);
    end
    % 生成非均匀相移的矩阵
    Alpha = cell(F,1);
    Delta_mac = Delta;Q = 1;
    for f = 1 : F
        Alpha{f,1} = [0, f, f, f, f, f] * phase_shift_tilts;
        [~,Delta_mac{f,1}] = Maclaurin(I_row,I_col,Q,Alpha{f,1});
        Delta_mac{1,1}(:) = 0;
    end
    
%     [~, Delta_lrd] = create_lrd(I_row,I_col,Q,phase_shift_tilts);

    % 均匀相移+非均匀相移 = 最终加在物体上的相移
    for f = 1 : F
        Delta{f,1} = Delta{f,1} + Delta_mac{f,1};
        Delta{1,1}(:) = 0;
        addDelta = Delta{f,1};
        save([filepath,'-',int2str(f),'-addDelta.mat'],'addDelta');
        mean_phase_shift = sum(addDelta(:))/(I_row*I_col);
        save([filepath,'-',int2str(f),'-mean_phase_shift.mat'],'mean_phase_shift');
    end

    
    % 生成物体
    d = 10;
    l = d*10;
    pk = peaks(I_row);
    pk = abs(pk);
    pk = 5*pk;
    pk = 2*pi*d*pk./((l-pk)*I_col/freq);

    for f = 1:F
        delta = Delta{f,1};
        for c = 1:I_col
              Phi = c*2*pi*freq/I_col.*ones(I_row,1) + pk(:,c);
              Figure{f,1}(:,c) = A + B*cos(Phi+mod((f-1),4)*pi/2+delta(:,c));
        end
    end

    num = 1;
    for f = 1 : F
       str = [filepath,'-',int2str(f),'.bmp'];
       Figure{f,1} = imnoise(Figure{f,1},'gaussian',noise_mean_num,noise_intensity);
       imwrite(Figure{f,1},str);fprintf('------%s条纹图生成！------ \n', str);
       num = num + 1;
    end
end